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Abstract: At the wall in a hybrid nematic cell with strong anchoring, the nematic director 
is parallel to one wall and perpendicular to the other. Within the Landau-de Gennes theory, 
we have investigated the dynamics of s = ±1/2 wedge disclinations in such a cell, using the 
two-dimensional finite-difference iterative method. Our results show that with the cell gap 
decreasing, the core of the defect explodes, and the biaxiality propagates inside the cell. 
At a critical value of d c * ~ 9b, (where t, is the characteristic length for order-parameter 
changes), the exchange solution is stable, while the defect core solution becomes 
metastable. Comparing to the case with no initial disclination, the value at which the 
exchange solution becomes stable increases relatively. At a critical separation of d c ~ 6^, 
the system undergoes a structural transition, and the defect core merges into a biaxial layer 
with large biaxiality. For weak anchoring boundary conditions, a similar structural 
transition takes place at a relative lower critical value. Because of the weakened frustration, 
the asymmetric boundary conditions repel the defect to the weak anchoring boundary and 
have a relatively lower critical value of d a , where the shape of the defect deforms. Further, 
the response time between two very close cell gaps is about tens of microseconds, and the 
response becomes slower as the defect explodes. 
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1. Introduction 

Topological defects arise as a result of broken continuous symmetry and are ubiquitous in nature, 
from microscopic condensed matter systems governed by quantum mechanics to a universe in which 
gravity plays a decisive role [1-3]. Defects in liquid crystals (LCs) have been the subject of much 
interest, still offering unsolved problems. Commonly observed defects in the uniaxial nematic phase 
are typically point defects with topological charge 5=1 and line defects with topological charge 
s = ±1/2 [4]. There are two types of s = ±1/2 disclination lines [5]. The first type is wedge disclination, 
with the rotation vector parallel to the disclination line, while the second is twist disclination, with the 
rotation vector perpendicular to the disclination line. The region where the presence of a defect causes 
apparent deviations from bulk ordering is referred to as the defect core [6]. 

Eigenvalue exchange (also called order reconstruction) was first shown by Schopohl and Sluckin 
within the core of s = ±1/2 wedge disclinations [7] and has aroused numerous subsequent studies on 
the detailed core structure in nematic LCs for various boundary conditions theoretically [8-13] and 
experimentally [14-17]. In particular, Palffy-Muhoray [8] investigated the eigenvalue exchange 
problem in a hybrid alignment nematic (HAN) cell with perpendicular alignment of NLC molecules on 
one surface and parallel alignment on the other. A structure transition from the eigenvector rotation 
configuration to the eigenvalue exchange is expected when the cell gap is decreased. Similar results 
are achieved in a twisted nematic (TN) cell with n/2 boundary twist configuration [9]. 

Recently, the atomic force microscope (AFM) has been used as a unique tool to measure the 
separation dependence of surface forces with unprecedented accuracy and flexibility [18-26]. This is 
not only important for studying the molecular and surface forces due to their potential application in 
self-assembly at the nanoscale, but is also important for the structural analysis of the interfaces of 
soft matter and complex fluids in contact with a solid wall. As an example, the structure of LC layers 
have been successfully analyzed using the force spectroscopy mode of an AFM [18-20]. Particularly, 
Carbone et al. [18] give the experimental evidence of the eigenvalue exchange in a thin HAN 
layer in the presence of a topological defect. We know that the speed of the AFM probe is about 
1-5 um/s, which is much slower, corresponding to the dynamic response of the structural transition 
(-10-100 us) [14-16], in other words, the speed of the AFM probe is slow enough for the system to 
undergo a structural transition; so, the AFM can be used to investigate the dynamic process of the 
structural transition as the call gap varies. 

Recently, many theoretical investigations have focused on the dynamical behavior of structure 
transition, using the Q tensor numerical model [27-31], which also describe the effect of order 
reconstruction-eigenvalue exchange. In this study, we will investigate the dynamics of biaxial 
transition in a nanoconfined HAN layer with a topological defect numerically, using the 
two-dimensional finite-difference iterative method. Our study is based on the Landau-de Gennes 
theory [32] describing the orientational order of an LC in terms of a second-rank tensor Q, which 
encompasses both uniaxial and biaxial states. 
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2. Theoretical Basis 

2.1. Free Energy 

Our theoretical argument is based on the Landau-de Gennes theory [32], in which the orientational 
order of an LC is described by a second-rank symmetric and traceless tensor [33]: 

e=iv/®«/- (i) 

1=1 

Here, X[ and e; are the z\h eigenvalue and the z' t h eigenvector of Q, respectively. In the isotropic phase, 
Q vanishes. In the uniaxial ordering, Q has two degenerate eigenvalues and can be represented by 

Q = S{?>n®n-\), (2) 

where S is the uniaxial scalar parameter, and the unit vector, n , is the nematic director pointing along 
the local uniaxial ordering direction. In Equation (2), 5* can have either sign: when it is positive, the 
ensemble of molecules represented by Q tends to be aligned along n , whereas when S is negative, 
it tends to lie in the plane orthogonal to n . 

Finally, when all eigenvalues of Q are distinct, the LC is in a biaxial state. The degree of biaxiality 
is measured by the biaxiality parameter, [3 , defined as [34] 

which is a convenient parameter for illustrating spatial inhomogeneities of Q and ranges in the interval 
[0, 1]. In all uniaxial states with two degenerate eigenvalues, p 2 = 0, while states with maximal 
biaxiality correspond to p 2 = 1. Since tr(Q 3 ) = 3detg, the states with [3 2 = 1 are precisely those where 
detQ = 0, which further implies that at least one eigenvalue of Q vanishes. 

Following the notation in [7], the Landau-de Gennes free energy density of LC is given by 
/=/buik{gap} +/e{g a p, v }, in which 



/h* = AtrQ 2 +^BtrQ 3 +^c{trQ 2 J (4) 

is the bulk energy that describes a homogeneous phase. In /bulk, B and C are positive constants, and A is 
assumed to vary with temperature Tin the form A = a {T - T*), where a is a positive constants and T* 
is the nematic supercooling temperature. Equation (4) gives the bulk equilibrium value of the uniaxial 

B I 24AC 

scalar order parameter in Equation (2), s eq = — — (1 + II — — ) , which also depends 
the temperature. 

The free-energy, f e , which penalizes gradients in the tensor order parameter field, is given in 
the form 

dx k dx k dxj dx k dx k dxj y ' 



Is on 



Int. J. Mol. Sci. 2013, 14 



24138 



The coefficients, L\, are related to the splay, twist and bend elastic constants. In this work, we use 
the one-elastic-constant approximation just for simplicity, where the splay, twist and bend elastic 
constants have a common value that depends quadratically on the scalar order parameter K = 2S 2 L U 
and Equation (5) reduces to 

The boundary conditions at the LC interfaces are taken into account by the surface term, f s , 
describing the interaction between the nematic molecules close to the substrate and the substrate itself, 
which is given by 

f s =^HQ-Qj, (7) 

where W s = W/S is the anchoring strength, with W being the anchoring strength in the Frank elastic 
theory, and Q s is the value of the tensor order parameter preferred by the surface [27]. 
The stable solution of Q on the boundary condition must satisfy the Neumann condition 

d/ s ... ¥ _ n 

where Vk is the &th component of the outward vector normal to the substrate. This expression allows us 
to impose on the surface substrate uniaxial, as well as biaxial conditions. For rigid anchoring, this is 
equivalent to imposing the Dirichlet condition on the surface, i.e., 

Q = Q S - (9) 

After an appropriate rescaling of the variables, we can simplify the calculating process. Here, 
we follow the rescaling of Schopohl and Sluckin [7] by defining the following dimensionless 

quantities: / = //[^ 4 /(9C)3] , Qij = ^ lj /q 0 > d = where q 0 = — ^ is the order parameter at 

the isotropic-nematic phase transition point and ^ = = ls the characteristic length for 

order-parameter changes; then, the rescaled form of Equations (4)-(7) can be written as 



3 ~ (10) 

7,=f^(e-eJ, <i2) 

where the rescaled parameter, A = defines the temperature scale and W s = with W 0 = —q 0 B%. 
We notice that the isotropic-nematic transition takes place at A 0 — -. 

In the scaling employed above, the rescaled uniaxial ordering has the form 
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Q = S(?>n®n-\), 



(13) 



where 5 = ^ = ^(1 + ^jl — ^) is the rescaled uniaxial scalar parameter at equilibrium. 



2.2. Geometry of the Problem 

We choose a HAN cell with strong anchoring boundary conditions; the nematic director is parallel 
to one wall and perpendicular to the other (see Figure 1). The model has two stable degenerate 
configurations with a rotation difference of n. The plates are placed at z = ±d/2 of a 
Cartesian coordinate system. The lengths, d x and d y , of the cell along the x- and the j-axes are much 
larger than d (d x ~d y » d). 

We study the structure of s = ±1/2 wedge disclinations. The disclination line is parallel to the y-axis 
of the simulation cells, and we seek solutions independent of y. 

At the two plates z = ±d/2, we enforce the uniaxial anchoring represented by (p-dn = 0 and 
<p+d/2 = 7t/2, that is 



Q- d/2 = s 
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(15) 



On the lateral walls at x = ±d x /2, we also prescribe fixed boundary conditions with uniaxial 
ordering. The total rotation at x = d x /2 is -jt/2, while at x = -dJ2 is nil. These boundary conditions 
are compatible with the generation of line defects with topological charge s = ~\I2, with which we 
shall also apply to line defects with topological charge s = 1/2, with the conditions at the lateral 
walls exchanged. 

Figure 1. The geometry of the problem. 
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2.3. Numerical Methods 

We compute the evolution of the LC with a dynamic theory for its tensor order-parameter field, 
Q{r,t). The local values of the scalar order parameter, 5*, and the director, n , can be obtained from Q 
through its largest eigenvalue and its associated eigenvector, respectively. According to [17], the 
evolution equation describing the dynamics of Q can be written as 



dQ 
dt 



J - + -tr 



SQ 3 \8Q 



(16) 



The coefficient, F, is given by F = 6D*/[1 - 2>tr(Q 2 )f, where D* is the rotational diffusion for the 
nematic. In these equations, it is assumed that Sf/SQ has been symmetrized. The stable solution of 
Q must satisfy the system Equations (16) and, on the boundary, the Neumann condition (Equation (8)) 
or Dirichlet condition (Equation (9)). 

The numerical calculations are to be intended with respect to the scaled variables. When the 

functional derivatives in Equation (16) are evaluated, one obtains a partial differential equation for Q 
in the bulk, 



dt 

with the boundary conditions, 



^ = -fllAQ + 2Q 2 + hr(Q 2 )(Q + 3l)- 2V 2 g 



(17) 



V z g(- d/2) = [q{- d/2)- a (- d/2)\ (18) 

V Z Q(- d/2) = -%- [Q{- d/2) - Q s (d/2)\ (19) 

where T = T x (— Bq 0 ) and V z = — . In general, Equation (17), subject to boundaries Equations (14) 

and (15), has more than one solution. One is that the structure satisfies the initial conditions given 
above, i.e., the configuration with a line defect; the second one is a pure HAN configuration consisting 
of a rotation of the eigenvector of Q (the director) with the z-axis only, i.e., a HAN configuration with 
no defect; and the third is the configuration with a biaxial layer, where there is little or no net rotation 
of the eigenvector, i.e., eigenvalue exchange configuration. In this work, we focus on the equilibrium 
configuration and dynamics with the initial conditions given in Section 2.2. 

We use a two-dimensional finite-difference iterative method employed in our previous studies 
in [35]. To discretize Equation (17), we replace the derivatives with finite difference as follows: 

t + At)-Q{t)]/At = -f ^2AQ{t)+ 2Q 2 {t) + hrQ 2 (t)(Q{t)+ 3l)- 2V 2 ^)}- (20) 

In our numerical calculations, we have found that a discretization with a time step given by 
4 x 10~ 10 is sufficient to guarantee the stability of the numerical procedure. In addition, our 
equilibration runs take 10 6 , which has been confirmed as sufficient for the system to reach 
equilibrium state. 
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We let the system relax from an initial condition under the boundary conditions given in Section 2.2. 
We have calculated the tensor Q after the system to reach equilibrium state, and then, the current 
tensor is diagonalized. Accordingly, the director is identified by the eigenvector possessing the largest 
eigenvalue. Then, the degree of biaxiality is measured by the biaxiality parameter, |3 2 , given by 
Equation (3), and the corresponding free energy is calculated. When the equilibrium state with a 
different cell gap is reached, we compute the dynamic evolution process between different cell gaps 
with the evolution equation given by Equation (16). 

3. Results 

In this section, we present our numerical results. According to the parameters given in [36], we 
have a = 0.043 x 10" 6 J/m 3 , B = -1.06 J/m 3 , C = 0.87 J/m 3 , L x = 2.25 x 10" 12 J/m. In our simulations, 

we set the scaled temperature at A = ^, corresponding to S = ^(1 + .JV^)- The rotational diffusion, 

D*, is set to 0.35, equal to the value used in [17]. We varied the value of d x to identify the behavior of 
the equilibrium configuration in the limit of d x — > qo. Our simulations are based on the initial 
conditions given in Section 2.2. We first focus on the strong anchoring conditions on the bounding 
plates and then on the weak anchoring conditions. 

3.1. Strong Anchoring Boundary Conditions 

3.1.1. The Structure at Equilibrium State 

Figures 2 and 3 show the director field profile and the calculated biaxiality, |3 2 (Equation (3)), inside 
the cell for different values of cell gap d. We find that at relatively large separations, the core structure 
is weakly influenced by d (see Figure 2a,b), and the nematic is uniaxial everywhere (|3 2 = 0), except for 
a small region around the defect core (Figure 3a,b). Its size is a few characteristic lengths, When the 
gap becomes comparable to the core of the defect explodes obviously along the x direction 
(Figure 2c) and a large biaxiality propagates inside of the cell (Figure 3c). The further decrease of the 
cell gap results in the creation of a biaxial wall, connecting the two orthogonal uniaxial directions 
imposed by the boundary conditions (Figures 2d and 3d) at a critical value of d c ~ 6£. In a word, the 
system has transited from the eigenvector rotation configuration with a topological defect into the 
eigenvalue exchange configuration by developing a thin transitional biaxial nematic layer. 

A more detailed analysis of the defect core, directly showing the process of biaxial transition, is 
given in Figure 4. There, we plot |3 across the defect center along z = 0 for different values of cell 
gap d. As the confining gap decreases, the distance at which |3 2 reaches its maximum value p 2 = 1 (i.e., 
the biaxial region) increased, indicating that the defect core becomes increasingly prolate in the 
x direction. By comparing the two figures in Figure 4, we see that at the cell gap where the defect core 
explodes, the changes of the biaxial region is faster. 
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Figure 2. The director field profile at equilibrium state in a thin nematic layer with 
different thicknesses. The simulations are based on the initial conditions given in 
Section 2.2. (a) d = 15% (b) d = 10% (c) d = 7.5% and (d) d = 6% 
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Figure 3. Biaxiality |3 2 for different thicknesses in a thin hybrid alignment nematic (HAN) 
layer with an initial line defect, (a) d= 15%; (b) d = 10% (c) d = 7.5% and (d) d = 6% 
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Figure 4. Biaxiality |3 2 across the defect center along z = 0 with different thicknesses, 
(a) Relatively large cell gaps; and (b) the cell gaps where the defect core explodes. 
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To quantify the influence of the cell gap on the defect structure, we plot in Figure 5 the change in 
both the width, r x , and height, r z , of the defect core as the cell gap, d, changes. We formally define r x 
and r z as the distance along the x- and z-axes, respectively, between the defect center and the point 
where P = 1. It is shown that r z is only slightly affected, while r x increases rapidly as the defect 
explodes. At a critical value of about d ~ 6^, the value of r x diverges approximately, signaling a 
nematic order reconstruction. Therefore, the contours representing the biaxiality |3 2 become 
increasingly prolate in the x direction and, finally, into a plane (see Figure 3). 

It seems that the transition to the biaxial wall originates from the biaxiality seed concentrated at the 
-1/2 defect. The biaxial wall would also be created in the absence of a defect [8,9]. According to [8,9], 
the two orthogonal uniaxial directors prescribed on the plate (HAN cell or TN cell with a n/2 twist) can 
be connected either through a director rotation or through an order reconstruction. When the cell gap is 
greater than a critical value d > d 0 o, the rotation solution is stable, while the exchange solution is 
unstable. When d < d c o, only the exchange solution exists as a stable configuration; no rotation solution 
exist. In other words, the rotation solutions merge continuously into the exchange solution at a critical 
value, d c o. 
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Figure 5. Variation of the defect core size with different thicknesses. The lengths of r x and 
r z estimate the core size along the x- and z-axes, respectively. 
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In order to make the mechanism of the structure transition clear, we make a comparison between 
the two cases (with and without initial disclination). Figure 6 gives the free energy variation inside the 
cell for three structures (the structure with a defect core, eigenvalue exchange structure and HAN 
structure with no defect), where the free energy is calculated from numerical simulation results. Here, 
the HAN structure plays the role of a case of comparison. 



Figure 6. The free energy as a function of d/t, for three different structures, (b) is a partial 
enlargement of (a). 
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Figure 6 shows that when the gap decreases to a critical value of d* ~ 9^, the exchange solution is 
stable, while the defect core solution is metastable. Comparing to the case with no initial disclination, 
we find that for the configuration with a topological defect, the value at which the exchange solution 
becomes stable increases relatively. That is because the configuration with a topological defect has 
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larger energy, relative to the pure HAN configuration with no defect. With further decreasing of the 
cell gap until a critical value, d c , a transition of the defect core solution into the exchange solution 
happens. Below d c , only the exchange solution exists, but no defect core solution. According to the 
structure given by Figures 2 and 3, we conclude that for the case with initial disclinations, the value at 
which only the exchange solutions exist is about d c ~ 6t,. 

Figure 6b presents that at a value of d c ~ 7^, the two cases (with and without initial disclination) 
have nearly the same energy, while the simulated transition value with initial disclination is at d c ~ 6t,. 
As mentioned in Section 2.2, we have prescribed fixed boundary conditions on the lateral walls, which 
will have an impact on the structure. In addition, there may be some calculation error. We speculate 
that for the two cases (with and without initial disclination), the values at which the structural 
transition takes place (only the exchange solutions exist) are almost the same, i.e., d c = d c o. 

3.1.2. Dynamical Evolution 

In order to make it clear how the configuration changes between two states, we make dynamical 
calculations of the defect core structure between different cell gaps (10^ — > 9^, 9t, — > 8^, 8^ — > 7^, 
7^ — > 6£). Figure 7 shows the evolution process of the calculated biaxiality |3 2 across the defect center 
along z = 0 as the cell gap changes. It is shown that if the call gap, d, is greater than the critical value 
of d c , an initial structure with a defect core cannot change to the eigenvalue exchange structure, even at 
the value where the exchange structure is stable. The response time is about tens of microseconds, and 
the response becomes slower as the defect explodes (see Figure 7c). 

To analyze the influence of the defect core on the dynamic process, we make dynamical 
calculations of the exchange configurations between different cell gaps. To aid the comparisons, we 
give the dynamics as the cell gap decreases. Figure 8 shows the evolution process of the calculated 
biaxiality, |3 2 , along x = 0 as the cell gap changes. It is shown that the response time is about several 
microseconds, and the response becomes faster as the cell gap decreases. Comparing Figure 7 with 
Figure 8, we find that the response between two exchange structures is much faster than that between 
two defect core structures. 

Figure 7. The dynamics of the biaxiality, p , across the defect center along z = 0 for the 
structure with a defect, (a) 10^ %; (b) % 8^; (c) 8^ 7^; and (d) 7^ 6^. 
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Figure 7. Cont. 




Figures 7 and 8 show that above the critical value, d c , there are only defect-defect and exchange- 
exchange transitions, but no defect-exchange or exchange-defect transitions existing. However, the 
response of an initial defect configuration with d > d c to d < d c (e.g., 7^ — » 6Q will undergo a 
defect-exchange transition (See Figure 7d). We expect that there are cell gap induced hysteresis effects, 
which are similar to hysteresis effects induced by the voltage in blue -phase liquid crystal [37]. As the 
cell gap, d, decreases, the defect core becomes increasingly prolate and undergoes a structural 
transition to the exchange configuration at a critical value of d c ~ 6t,. While the cell gap, d, is 
increasing from d < 6^, the system is in the exchange configuration until d ~ 9^, above which the 
system will change to the defect core configuration under the effect of thermal disturbance. 

Figure 8. The dynamics of the biaxiality, P , along x = 0 for the exchange solution, 
(a) m - 9& (b) % - 8^; (c) 8^ - 7^; and (d) 7^ - 6^. 
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3.2. Weak Anchoring Boundary Conditions 

In order to get the influence of the surface anchoring at the boundary plates on the behavior of 
topological defects, we focus on the Neumann conditions in this section. First, we study the condition 
of asymmetric boundaries. We impose strong homeotropic anchoring on the top plate and weak planar 
anchoring on the bottom plate. Figure 9 shows the director field profile inside a cell of thickness 
d = \0t, for different values of W s at the bottom plate. It is shown that as the anchoring strength, W s , 
decreases, the defect glides towards the weak anchoring bottom plate and finally coalesces with the 
surface layer. This is consistent with the results in [38,39]. 

Figure 9. The director field profile at equilibrium state in a thin nematic layer with 
different anchoring strength on the bottom plate. The cell gap of the simulated nematic 
layer is d= 10^. (a) W s =3;(b) W, = U(c) W, =0.5; and (d) W s =0S. 
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Figure 9. Cont. 
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Next, we fix the anchoring strength at the bottom plate to W s = 1, to investigate the influence of the 
cell gap on the topological defects. Figures 10 and 11 show the director field profile, as well as the 
change in both the width r x and height r z of the defect core inside the cell for different values of 
cell gap d. 

Figure 10 shows that the defect core is not located in the middle of the cell, which is due to the 
asymmetric boundary conditions. As the cell gap, d, decreases, the line defect moves towards the 
bottom plane. At a critical value d c (w as ) ~ 5£, the defect is transformed into a surface layer. The 
surface plane defect is then expelled from the cell when the cell gap reaches d c **(w as ) ~ 2.5£,. 
For d c **(w afi ) < d < d c (w as ), the surface biaxial layer bridges two uniaxial stares, one homeotropic in 
the bulk and the other planer on the boundary. The results are consistent with the results in [12,27]. 

Figure 10. The director field profile at equilibrium state in a thin nematic layer with 
different thicknesses. The anchoring strength is W s = 1 at bottom plates, (a) d = 10^; 
(b) d = 7^; (c) d = 5^; (d) d = 4^; and (e) d = 2.5£ 
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Figure 11. Variation of the defect core size with different thicknesses for the asymmetric 
boundaries. The lengths of r x and r z estimate the core size along the x- and z-axes, 
respectively. In the inset, the small region is enlarged. Here, the strong anchoring boundary 
plays the role of a case of comparison. 
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Figure 1 1 shows that the change in r z is slight, while r x increases when the cell gap reaches a critical 
value, d a , where the defect core begins to lose its circular cross-section and stretches along the 
x direction. By comparing with the strong anchoring conditions (see the inset in Figure 11), we can 
conclude that the asymmetric boundary condition has a relatively smaller value of d a . That is because 
the moving of the defect as the cell gap decreases weakens the frustration on the confining surfaces. 
Figure 11 also gives that the value, d c (w as ), where the defect is transformed into a surface layer, is 
smaller than that of strong anchoring conditions. 

Further, we study the condition of symmetric boundaries. We impose the same value of weak 
anchoring on both the two boundary plates and studied the effect of anchoring strength on the critical 
value, d c . In our simulation, the anchoring strength is set equal to W s = 1 at both plates. Figure 12 
shows the director field profile inside the cell for different values of cell gap d, which shows that the 
structural transition takes place at a critical value of d c (w s ) = 3£,. Compared with the strong anchoring 
conditions (i.e., Figure 2), we can conclude that weak anchoring boundary conditions decrease the 
critical value, d c . That is because the weak anchoring boundary weakens the frustration on the 
confining surfaces. 

Figure 12. The director field profile at equilibrium state in a thin nematic layer with 
different thicknesses. The anchoring strength is W s = 1 at both plates, (a) d = 6^; 
(b) d = 4!~; (c) d = 3.5^; and (d) d = 3£ 
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4. Conclusions 

Within the Landau-de Gennes theory, we carried out a numerical study on the structure of s = -1/2 
wedge disclinations and investigated the dynamics of biaxial transition in a nanoconfined HAN layer, 
using the two-dimensional finite-difference iterative method. Our results show that with the cell gap 
decreasing, the core of defect explodes and the biaxiality propagates inside the cell. Below a critical 
thickness, d c , of a HAN layer with a topological defect, an eigenvalue exchange transition occurs, and 
the defect is squashed flat into a biaxial layer. Because of the presence of the defect core, the value at 
which the exchange solution becomes stable increases relative to the pure HAN. The dynamic 
simulations show that before the critical value of d c , an initial structure with a defect core cannot 
change to the eigenvalue exchange structure, even at the value where the exchange structure is stable. 
The response time between two very close cell gaps is about tens of microseconds, and the response 
becomes slower as the defect explodes. For weak anchoring conditions, plates with strong anchoring 
repel the defects, while plates with weak anchoring allow the defect to escape through the boundary. 
Because of the weakened frustration, the weak anchoring boundary condition has relatively lower 
critical values of d a and d c , corresponding to the value where the defect core begins to lose its circular 
cross-section and where the defect is transformed into a surface layer, respectively. 

According to the parameters given, we can easily get the characteristic length ^ ~ 3.96 nm, and the 
critical separation is -20 nm, which is consistent with the experimental result in [18]. Furthermore, the 
response time also agree with the experimental results (-10-100 us) in [14]. The consistency of 
numerical calculation and experiment confirms the correctness and rationality of our theory, which 
further illustrates that the AFM can be used to carry on dynamics research. 

Further, the parameters also give that the unit of anchoring strength is Wo ~ 5.7 x 10~ 4 J/m 2 , which 
is reasonable. According to b = k/W, W s = 0.1 corresponds to an extrapolation length of b ~ 10~ 7 m, which 

is larger than the cell gap of d = \0t, ~ 3.96 x 10~ 8 m, i.e., bid > 1. Figure 9 shows that at W s =0.1 , 

a sample of d = 10^ presents a homeotropic orientation, which is consistent with the results of [39]. 

It can be predicted that the critical separation and the response time should be affected by the temperature; 
the detailed results, as well as the force curve with the Q tensor numerical model are tasks for the future. 
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